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Abstract 

We obtain iso-spectral Euler-Bernoulli beams by using factorization and Lie sym- 
metry techniques. The canonical Euler-Bernoulli beam operator is factorized as the 
product of a second-order linear differential operator and its adjoint. The factors are 
then reversed to obtain iso-spectral beams. The factorization is possible provided 
the coefficients of the factors satisfy a system of non-linear ordinary differential 
equations . The uncoupling of this system yields a single non-linear third-order 
ordinary differential equation. This ordinary differential equation, refer to as the 
principal equation, is analyzed and solved using Lie group methods. We show that 
the principal equation may admit a one-dimensional or three-dimensional symmetry 
Lie algebra. When the principal system admits a unique symmetry, the best we can 
do is to depress its order by one. We obtain a one-parameter family of solutions 
in this case. The maximally symmetric case is shown to be isomorphic to a Chazy 
equation which is solved in closed form to derive the general solution of the principal 
equation. 



1 Introduction 



A recent study [1] suggests that efforts to model the transverse motion of vi- 
brating beams date back to Leonardo da Vinci. In his discussion of the bending 
of beam/ spring with rectangular cross-section he wrote [2]: "Of bending of the 
springs: if a straight spring is bent, it is necessary that its convex part become 
thinner and its concave part, thicker. This modification is pyramidal, and con- 
sequently, there will never be a change in the middle of the spring". About 
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hundred years after Da Vinci's attempt to develop a beam theory, Galileo 
suggested an erroneous calculation of the load carrying capacity of a trans- 
versely loaded beam. Da Vinci did not benefit from Hooke's law and calculus 
which postdate him whereas Galileo made the incorrect assumption that under 
transverse loading, the stress is uniformly distributed on cross-sections. The 
first correct and systematic formulation of an elasticity theory was done by 
Jacob Bernoulli (1700-1782). He found that the curvature of an elastic beam 
at a given point is proportional to the bending moment at that same point. 
Taking advantage of his uncle seminal work, Daniel Bernoulli (170 7- 1783) de- 
rived the partial differential equation that governs the motion of a vibrating 
beam. Leonard Euler(1707-1705) pursued the study of the Benoullis by using 
and extending their theory in his investigation of the shape of elastic beams 
subjected to various external forces. For a recent review of beam theories, the 
reader is referred to the paper by Han, Benaroya and Wei [3] and references 
therein. 

In Euler-Bernoulli beam theory, the transverse motion of a thin non-uniform 
beam is governed by the partial differential equation 



dx 2 



E(x) J(s)g 



2 u 

+ p{x) a(x)— — = 0, < x < L, 



where E is the modulus of elasticity (a.k.a Young's modulus), a(x) is the 
cross-sectional area at x, I(x) is the area moment of inertia of the cross-section 
located at x about a neutral axis, p is the density, and L is the length of the 
beam. Equation (1) is solved subject to appropriate boundary conditions at 
x = and x = L. Some common boundary conditions and their physical 
meaning are 



— — - = 0, u — (hinged end) , (2) 
ox 

du 

— — = 0, u = (clamped end), (3) 
d 2 u d / . w/ . d 2 u\ , . 

— = 0,-^)1(^=0 (free end), (4) 
du I d 2 u\ 

_ = 0,^(E W /M^)=0(s M in g end). (5) 

The ansatz u(t,x) = e iuJt Y(x), where u is the frequency of vibrations, solves 
Eq. (1) provided 

& 2 { f{x)d ^) =u2m{x)Y > °^ L > ( 6 ) 
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where 



f(x) = E(x) I(x), m(x) = p(x) a(x). 



(7) 



Under Barcilon's transformation [4,5] 



.t 



z = 



J {m/ffl\ Y = (m 3 /)- 1/S U, 



(8) 



o 



Equation (8) becomes 



d A U 



dU_ 

dz 



+ B(z)U = lu 2 U, 0<z<l, 



(9) 



dz 2 



where the correct expression of A and B are given by Gottlieb [5] . Save for the 
clamped end conditions, Barcilon's transformation leaves the other boundary 
conditions invariant under some constraints given explicitly in [5] . As we shall 
see in the sequel, the canonical form (9) plays a crucial in inverse spectral 
problems associated with Euler-Bernoulli beams. 

In many practical problems such as the structural identification through non- 
destructive means, one is interested in determining the physical characteristics 
of the beam i.e. f(x) and m(x) from few frequencies (spectra). That is, one 
tries to 'hear' the shape of the beam. Since the vibration of uniform beams (i.e. 
/ and m are constant functions) is well understood, an interesting problem 
related to the inverse spectral problem is to find non-unform beams that have 
the same length and share the same spectra with the uniform beam under the 
same boundary conditions. This problem was comprehensively studied by Got- 
tlieb [5] using the canonical form (9). He found seven classes of clamped beams 
that are iso-spectral with the standard unit-coefficient beam (A = = B). 
Some of these beams where rediscovered by Abrate [6] (see [7] for comments) 
using a different paradigm. For discrete beams, Gladwell [8], introduced two 
procedures for finding iso-spectral families. His first procedure relies on shifted 
orthogonal matrix triangulation whereas the second utilizes Toda flow. 

This paper is dedicated to the search of iso-spectral beams using Lie sym- 
metry methods together with a factorization method initiated by Pochel and 
Trubowitz [9] for Sturm-Liouville operator. To be more precise, the factoriza- 
tion method for the canonical beam operator leads to a system of nonlinear 
ordinary differential equation that we shall solve by Lie's reduction method. 
We have organized the present paper as follows. There are four section of 
which this introduction is the first. In section 2, we briefly motivate the use of 
factorization in the search of iso-spectral beams. The factorization approach 
leads to a system nonlinear ordinary differential equations. We employ Lie's 
symmetry theory to study its continuous symmetries in Section 3. Section 4 
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employs the symmetries calculated in Section 3 to find exact solutions of the 
system of nonlinear ordinary differential equations derived in Section 2. This 
results in new families of iso-spectral beams. 




2 Iso-spectral deformation through factorization 

The main idea is to factorize the canonical beam operator as the product 
of a linear second-order differential operator and its adjoint, then swap the 
order of the operators to obtain iso-spectral operators. This idea was already 
used for vibrating rods by Pochel and Trubowitz [9]. Ghanbari [10] attempted 
the factorization method on vibrating beams but failed to solve the resulting 
system of nonlinear ordinary differential equations. 

Let us define the beam operator as the differential operator 

+ B(z). (10) 

Factorize the differential operator £ as 

£ = ?m, (ii) 

where 

K =^ + r W^ + s W- *" = &-TJ r ^- ] + s ^ (12) 

r and s are smooth functions of their argument, and JH* is the adjoint of 
Simple reckoning shows that the factorization given in Eq. (11) constraints 
the functions r and s to satisfy the following system of nonlinear ordinary 
differential equations [10]. 

r' - r 2 + 2s = A(z) (13) 
s" -(rs)' + s 2 = B(z), (14) 

where the prime stands for differentiation with respect to z. 

Reversing the operators in the factorization (11), we obtained the operator 
[10] 



where 
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A(z) 

B(z) 



= 2s-3r'-r 2 , 

= s 2 + s" - r'" - rr" + rs' - sr'. 



(16) 
(17) 



Define an eigenpair of £ as being a pair (u,U), where U is not the zero 
function, such that 2XJ = uo 2 U . It can be easily verified that if (u, U) (resp. 
(u,U)) is an eigenpair of £ (resp. £), then (u,fRU) (resp. (lj,DVU) is an 
eigenpair of £ (resp. £). Thus, loosely (we have not considered the boundary 
conditions), £ and £ are iso-spectral. This motivates the use of factorization 
in the search for iso-spectral beam. 

In order to characterize iso-spectral beams obtained via factorization, one 
needs to solve Eqs (13)-(14) for r and s. Solving Eq. (13) for s, and substituting 
into Eq. (14) yield the single nonlinear ordinary differential equation for r 

7 A 2 r 4 
r (3) _ 3rr " _ _ r >2 + 2 ( 2r 2 + A y _ A "_ r 2 A + rA t + 2 B = 0.(18) 

In the next two sections, we shall use Lie symmetry analysis to solve Eq. (18). 
We shall refer to Eq. (18) as the principal equation. 



3 Lie symmetry analysis of the principal equation 



A symmetry of a differential equation is an invertible transformation of the 
dependents and dependents variables that leaves the equation unchanged. The 
main attractive property of symmetries is that they transform solutions into 
solutions, and better, they provide a systematic route to integration. 

Historically, the Norwegian mathematician Sophus Lie (1842-1899) realized 
that determining all the symmetries of a given equation is a formidable task. 
However, if we restrict ourself to symmetries that depends continuously on 
a small parameter and that forms a group (called an infinitesimal group in 
Lie's terminology or Lie group in modern terminology) , it is possible to use 
the machinery of calculus to determine such symmetries in an algorithmic 
fashion. Strikingly, Sophus Lie discovered that symmetries are road maps to 
integration. 

Here, we are concerned with the determination of the Lie symmetries of Eq. 
(18). For more details about Lie's symmetry theory, the reader is referred to 
the books [11,12,13]. 

A vector field 

X = ^ z i r )ir + r i{^r)—, 



(19) 



is a point symmetry of Eq. (18) if 
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r (3) _ 3rr " _ Z r '2 + 2(2r 2 + _ A „ _ r 2 A 



A 2 r A 



= 0. (20) 

Eq.(i8) 



In Eq. (20), the notation |gq ^ lg ^ means that after expanding the left hand 

side, Eq. (18) has to be used to eliminate r^ 3 \ The operator is the third 
prolongation of X defined recursively by 

X [k+i] =x [k] + [k]d^ X [ 0]=X) (21) 

Of ) 

v [k] = D ^[fc-i]) + r WD(£), r] [0] = 77, (22) 
i9 Q 

dz dr dr' ^ ^ 

Since the symmetry coefficients £ and rj are independent r' and r", Eq. (20) 
is polynomial in these derivatives. Thus we may set coefficients of the mono- 
mials (r') m (r") n to zero in Eq. (20). It results an over-determined system of 
linear partial differential equations for £ and rj that simplifies to the following 
equations. 

£ = a(z), r ] = -a'(z)r-2a"(z), (24) 

aA' + 2d A + 5a (3) = 0, (25) 

r,/ 4a ' r, «' ,2 o!', 2a' „ 2a< 3 > A® a< 5 > 

£' + — -B = — A H — - — | — j A + — A + A H — - — h 26 

aa 2 a' a a 2a 

3a' A' + 2a"A + aA" + 5a (4) = 0, (27) 
where a is an arbitrary smooth function of z. Simple calculations show that 

^( Eq.(25)) = Eq.(27). 

Thus Eq. (27) is a mere differential consequence of Eq.(25). In order to solve 
the remaining equations, we consider the following cases. 

Case I : A = and B = 0. In this case, we find that a = k\Z 2 + k^z + k 3 , 
where ki, k 2 , and k 3 are arbitrary constants. The symmetry Lie algebra is 
spanned by the operators 

*i = 7T» X2 = z TT- r TT^ X 3 = ^ 2 |-- 2(^ + 2)|- • (28) 
oz oz or oz or 

The Lie bracket of these operators are: 

[Xl, X2] = Xi, [Xi,X 3 ] = 2X2, IX2, X3] = X3. 
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The symmetry Lie algebra is non-solvable since its derived algebra of any order 
is nontrivial. To be more precise, if & =< Xi, X 2 , X 3 >, then &^ = &. In 
fact, 6 sZ(2,R) [14]. 

Case II: A ^ or B ^ 0. Solving Eqs. (25)-(26), we obtain after some 
calculations 

A= 5 "' 2 - 1 °°f + 2C ', (29) 
81a' 4 + 12a ,2 (3d - 17aa") + 72a 2 a'a< 3 ) 

£5 = 



16a 



4 



4 (C? + AC 2 - 6C ia a" + 21a 2 a" 2 - 6a 2 a< 4 >) 



+ ^ m? L • (30) 

where C\ and C 2 are arbitrary integration constants. The symmetry Lie alge- 
bra is spanned by the single operator 

r = a — - (a'r + 2a") — ■ (31) 
oz or 



4 Solutions of the principal equation and iso-spectral beams 

Here we exploit the symmetry structure of the principal equation (18) to find 
its solutions. 



4.1 A and B are given by Eqs. (29)-(30) 

In this case the symmetry Lie algebra is one-dimensional. The best we can 
do using Lie's integration technique is to reduce the order of the equation 
by one. This is accomplished by introducing new dependent and independent 
variables as independent solutions (known as basis of first order differential 
invariants ) of 

r [1] / = 0. (32) 

A fundamental set of solutions of Eq. (32) are obtained by solving the system 
of ordinary differential equations 

dz dr dr> (33) 



a(z) -a'r - 2a" -2aV - a"r - 2a( 3 ) 

The solutions of Eq. (33) are 

ar + 2a! = q u a 2 r' + aa'r + 2aa" = q 2 , (34) 
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where q\ and q 2 are arbitrary constants of integration. Thus we introduce the 
new variables 

u = ar + 2a', v = a 2 r' + aa'r + 2aa" . (35) 
In terms of the new variables, Eq. (18) reads 

2v 2 = 6uv ^l - 2v (^X + u 4 - 8u 2 v + 7v 2 + 2C lU 2 - AC x v - AC 2 .(36) 
au l du \du J 

Equation (36) is devoid of Lie symmetries so that we may not further integrate 
Eq. (36) by Lie's method. Assume from now on that CI = = C 2 , and look 
for a solution of Eq. (36) in the form 

v = ku 2 , (37) 

where k is a constant to be determined. Substituting the ansatz (37) into Eq. 
(36) yields 

12k 3 - 19k 2 + 8k - 1 = 0. (38) 
Solving Eq. (38), we obtain 

fc ^ (ill). (39) 



.4'3 

Equation (37) may be written as 

d ( a'\ a! ( a' \ , / a' \ 2 . . 

— \r + 2-)+-[r + 2-)=k[r + 2-). 40 
dz \ a J a \ a J \ a J 

Whence the natural change of variable 

w = r + 2- ■ (41) 
a 

In the new variable, Eq. (40) becomes Bernoulli's equation 
a' 

w ' + — w = kw 2 . (42) 



a 



The classical technique for integrating Bernoulli's equation leads to 

W = C a - ka j z a-^dt ' (43) 
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where C is an arbitrary constant of integration. Hence 
1 2a' , fl 1 I 

r = 7T—^-7r-n:-—, ke -,-,l . (44) 



Ca-kaJ'a^dt a ' U 3 

To summarize, we have established the following result. 
Theorem 1 TTie operator £, where 

5 a' 2 — 10 aa" , 

A = 2a2 ' (45) 
„ 81 a' 4 - 204 aa'V + 72 aVa® + 84 a 2 a" 2 - 24 aV 4 ) 

B = 16^ ' (46) 

iso- spectral with the operator SL, where A and B are defined by Eq. (16)- (17), 
r is given by Eq. (44) an d s — (A + r 2 — r')/2. 

^.^ Case of the standard unit- coefficient beam (A = and B = 0) 



According to the symmetry analysis performed in section 3, this is the most 
symmetric instance of the principal equation. It turns out that, in this par- 
ticular case, the principal equation can be invertibly mapped to a Chazy 
equation[15,16,17]. Indeed, Chazy's equation 

Vxxx = 2yy xx - 3yl + a(6y x - y 2 ) 2 , a = const., (47) 
admits a 3D symmetry Lie algebra spanned by the operators [18] 

* = * = Y ^ x2 l- {2xv + e) Vy- (48) 

Simple calculations show that < X 1: X 2 ,X 3 > and < Yi,Y 2 ,Y 3 > are two 
equivalent (up to an invertible transformation) representations of s/(2,R). 
Indeed the transformation 

2 . . 

r = y, z= z x i ( 49 ) 

maps < Xi,X 2 ,X 3 > into < Y 1 ,Y 2 ,Y 3 >. The same transformation maps the 
principal equation into Chazy's equation 

4 

y X xx = 2yy xx - 3y 2 x + — (6y x -y 2 ) 2 . (50) 

The integrability of Eq. (47) was established by Chazy [17]. The justification 
of the integrability of Eq. (47) from a Lie symmetry standpoint is due to 
Clarkson and Olver [18]. 
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The integration of Eq. (47) relies on the following theorem. 

Theorem 2 (Chazy [17]) Assume that if and ip are two arbitrary linearly 
independent solutions of the hypergeometric equation 

Then the general solution of Eq.(47) is given in parametric form by 

m' y ~mdx [ } 

Remark. The case a = 1/9 was treated by Clarkson and Olver [18]. In this 
degenerate case, the solution of Eq. (47) is still given by Eq. (52), where this 
time, if and ip are two arbitrary linearly independent solutions of the Airy 
equation x + ( c t/2)x — 0, an d c is an arbitrary constant. 

For Eq. (50), a = —1/48 and we may choose if and ip as [19] 

¥> = h 2 Fi (\, ~ §; 1 - t) - 4-^(1 - tf* 2Fl (^4 1 " *) 53) 

1> = h 2F 1 (1 -1 |;l-t) - 4-^(1 - t) 1/s 2F1 Q, ^ | 1 " ^ 

where 2-F1 is Gauss hypergeometric function, k\ to £4 are arbitrary constants, 
and the factor — 4~ 2 / 3 is introduced for convenience. Since the formula (52) 
is invariant under equal scaling of if and ip, we may assume without loss of 
generality that k\k 4 — k2k 3 = —1, so that Eqs. (53)-(54) involve only three 
arbitrary constants. Thanks to the exact formulas [20] 

(57) 

it is natural to introduce the re-parametrization 

1 4(2r-l)3 (58) 

Thus f and ip reduce to 

v = (k 1 + k 2 r l ^)(l-2r)- 1 /\ (59) 
i l J = (k 3 + k 4 r 1 ^)(l-2r)- 1 /\ (60) 
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Therefore the general solution of Eq. (50) is given in parametric form by the 
equations 

_ fcl + k 2 T^ _ 3{k, + k^){3hT^ + k A {2-T)) 

X ~h + k A TW V ~ l-2r ' 1 J 

Solving Eq. (61a) for r and substituting into Eq. (61b) yield 

_ 3 [3k 3 (k 3 x - ki) 2 (k 4 x - k 2 ) + 2k 4 (k 4 x - k 2 ) 3 + k 4 (k 3 x - kj) 3 ] 
V ~ (k 4 x - k 2 ) [2(kl - hxf - (k 4 x - k 2 f] ' { } 

Finally we obtain 

_ 6 [3fc 3 (3fc 3 ^ - 2k l ) 2 (3k 4 z - 2k 2 ) + 2k 4 (3k 4 z - 2k 2 f + k 4 (3k 3 z - 2h) 3 L 
T ~ {3k 4 z - 2k 2 ) [2(2fcl - 3A; 3 ^) 3 - {3k 4 z - 2k 2 f] ^ 

s =(r 2 -r')/2, (64) 
kik 4 — k 2 kz = —1. (65) 

Thus we have established the following result. 

Theorem 3 The unit beam operator -— is iso-spectral with the operator £ 

Uj/C 

with 



A = -5r', (66) 
B = -^(6r^ + 2rr" - r 4 - 7r' 2 ), (67) 

where r is defined by Eq. (63). 



5 Conclusion 



The factorization of the canonical beam operator as the product of a second- 
order linear differential operator and its adjoint results in a system of non- 
linear ordinary differential equations. We uncoupled the system of ordinary 
differential equations and solved the resulting nonlinear ordinary differential 
equation (viz. the principal equation) using Lie's method. The symmetry anal- 
ysis reveals that the principal equation admits either a one-dimensional or 
three-dimensional symmetry Lie algebra. When the principal equation admits 
a one-dimensional Lie algebra, its order can be reduced by one at most. In 
this case we obtained a one-parameter family of solutions (see Eq. (44)). For 
the maximally symmetric case, we proved that the principal equation can 
be mapped to a Chazy equation. The latter is solved in closed form and the 
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general solution (see Eq. (63)) of the principal equation is obtained. By revers- 
ing the order of the factorization of the canonical beam operator, numerous 
nontrivial iso-spectral families are obtained. 

It is opportune to mention that Nucci [21] calculated the Lie point symmetries 
of the system (13)-(14) without showing how these symmetries are used. The 
truth is, for systems of ordinary differential equations, successive reduction of 
order using symmetries can be ambiguous: there are in general several possi- 
bilities for choosing the new variables as invariants [22,23]. We avoided this 
problem by first uncoupling the system to obtain a single nonlinear ordinary 
differential equation on which we applied Lie symmetry analysis. 
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